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The successful prediction of light element abundances from Big Bang Nucleosynthesis (BBN) has 
been a pillar of the standard model of Cosmology. Because many of the relevant reaction rates are 
sensitive to the values of fundamental constants, such as the fine structure constant and the strong 
coupling constant, BBN is a useful tool to probe and to put constraints on possible cosmological 
variations of these constants, which arise naturally from many versions of extra-dimensional theories. 
In this paper, we study the dependences of fundamental constants on the radion field of the universal 
extra dimension model, and calculate the effects of such varying constants on BBN. We also discussed 
the possibility that the discrepancy between BBN and the Wilkinson Microwave Anisotropy Probe 
(WMAP) data on the baryon-to-photon ratio can be reduced if the volume of the extra dimensions 
was slightly larger - by 0(10"^) - at the BBN era compared to its present value, which would result 
in smaller gauge couplings at BBN by the same factor. 

PACS numbers: 



I. INTRODUCTION 



The standard Big Bang Nucleosynthesis (BBN) the- 
ory is a successful marriage between standard Friedmann 
cosmology and nuclear physics, explaining the origin and 
abundances of the light elements D, ''He, '^He and ^Li. 
Its only one input parameter, the baryon-to-photon ratio 
•q = nhjn-f (where rih and are respectively the num- 
ber density of baryons and photons) has now been deter- 
mined by the observations of the Wilkinson Microwave 
Anisotropy Probe (WMAP) with rather good accuracy 
0, so that there are essentially no free parameters in 
this scenario. Because the BBN predictions are sensitive 
to a set of physical quantities which depend on various 
fundamental constants such as the gauge couplings, the 
Yukawa couplings and gravitational constant (see Sec- 
tions III and IV for details), it would provide stringent 
constraints on the cosmolojgical variations of these con- 
stants. For example, in [3^13, Q the authors use BBN 
to constrain the change of fine structure constant, in 
constraints on the change of Higgs vacuum expec- 
tation value are considered, and the effects of a varying 
strong coupHng are discussed in e.^. H, I3i the varia- 
tion of gravitational constant has also been investigated 
[Ullllllllllll in the context of BBN. Besides, there is 
another reason for the interests in the interplay between 
varying fundamental constants and BBN: in contrast to 
the excellent consistency between BBN theory and ob- 
servation with T] given by WMAP for the deuteron (D) 
abundance, the predicted *He and ^Li abundances are 
smaller than the results implied by WMAP 15, 1(|; it is 
then suggested that the variations of fundamental con- 
stants such as the fine structure constant 0, ^3 or the 
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deuteron binding energy Bd 0| {Bd itself is certainly 
not a fundamental constant, but its variation possibly 
originates from some other fundamental constants as sug- 
gested in 0) eic. might partly or completely solve this 
discrepancy. 

The theoretical investigations of varying fundamental 
constants date back to the early work of Dirac in 1930's, 
and have re-aroused great attentions because of the re- 
cent discovery of Webb et al. |0| that the quasar 
absorption lines at redshifts of z = 1 ~ 3 suggest a small 
evolution of the fine structure constant between that pe- 
riod and present. Although variations of the fundamen- 
tal constants do not occur in the standard model (SM) of 
particle physics, the two leading paradigms of the physics 
beyond SM, namely the string theory '21] and the Kaluza 
Klein theory generally predict such variations. In 
these scenarios a dilaton field or the size of extra dimen- 
sions may evolve cosmologically; thus until their vacuum 
expectation values (VEV) are fixed, there might be a 
co-variation of several fundamental constants. This pic- 
ture is similar to the discussions in grand unified theories 
(GUTs) 01 and implies that considerations of one 
varying constant alone may be incomplete, although in- 
deed such a CO- variation itself is in general rather model- 
dependent (see [2^ for one case derived from string 
theory) . For a detailed discussion of the theoretical and 
observational (experimental) aspects of varying funda- 
mental constants, see p^ . 

The recent interests in the extra dimensional theo- 
ries have stimulated other considerations of varying fun- 
damental constants. Of the frequently discussed extra 
dimension models, the orig inal ADD 2^ and RS "2^ 
brane models do not induce changes in the gauge cou- 
plings even though their moduli fields evolve [S^, |3ll be- 
cause of the conformal invariance of the gauge kinetic 
terms (see however |32j for an alternative), and so in this 
work we will concentrate on the universal extra dimen- 
sion (UED) model, in which the fields can propagate in 
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all dimensions. This model was first proposed in |33| 
and became extremely interesting for cosmology since it 
was found later that it predicted the presence of a stable 
massive particle, the lightest KK partner (LKP), which 
is a natural dark matter candidate (see 34, and ref- 
erences therein for details). Our purpose in this work is 
twofold: firstly, we study how the fundamental constants 
in the low energy effective theory change if the size of 
the universal extra dimensions, the radion field, under- 
goes a slow cosmological evolution and to what extent the 
changes are allowed by the BBN observations; secondly, 
we show that if the volume of the extra space at the time 
of BBN is slightly larger than its present value, the dis- 
crepancy between BBN and WMAP discussed above may 
be reduced. 

The arrangement of this work is the following: in Sec- 
tion II we derive the effective low energy actions for the 
gravitational and matter sectors in the model; then we 
shall find the radion dependences of the fundamental 
constants in Section III. In Section IV we briefly review 
the standard BBN, point out the most relevant physical 
quantities for its prediction and relate these quantities 
to the radion field. Section V contains the numerical re- 
sults of this work, which are obtained by modifying the 
standard code of BBN [^Bl and including the effects dis- 
cussed in Section IV. Finally Section VI is devoted to 
discussion and conclusion. Throughout this work we as- 
sume 3 species of massless neutrinos as in the standard 
BBN and adopt the units h ~ c ~ 1 . 



II. THE LOW ENERGY 4-DIMENSIONAL 
EFFECTIVE ACTIONS 

Let us consider a general A+n dimensional model, with 
n being the number of the (universal) extra dimensions 
(although there may also be large extra dimensions, we 
do not consider them in this work) . The full line element 
is given as: 

ds^ = GABdX^dX^ = g^^dxf'dx'' + habdy''dy\ (1) 

where ii, v — 0, 1, 2, 3 label the four ordinary dimensions, 
a,6 = 4,---,3-|-n denote extra dimensions and A, B — 
0,l,2,---,3-|-ri the whole spacetime. For simplicity we 
shall not consider cross terms such as Gafi in Eq. (1). 
The extra dimensions are assumed to compactify on an 
orbifold, and their coordinates ya take values in the range 
[0, 1]. The quantities hab have dimensions of [Length]^ 
since ya are dimensionless in our choice. 

Because the energy range we are interested in is much 
lower than the inverse size of the universal extra dimen- 
sions, which is thought to be larger than several hundred 
GeV's, it is adequate to consider only the zero modes of 
the metric. Then the effective 4-dimensional action (in 
the gravitational sector) can be obtained by dimension- 
ally reducing Eq. (1) as: 
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(2) 



in which \g\, \h\ and \G\ are respectively the determinants 
of the metrics of the ordinary dimensions, the extra di- 
mensions and the whole spacetime. R^lg] and i?4+„[G] 
are the Ricci scalars of the ordinary 4 and the total A + n 
dimensional spacetimes. K4,K4+n are related to the 4 
and A + n dimensional Planck masses via = 2M| and 
= 2M^^^, while they themselves are connected by 
a volume suppression k|_|_„ = k| • with V being a 
measure of the extra space volume whose present-day 
value is denoted by Vq in Eq. (2) (Note that because 
of the specified choice of Vq and because the higher di- 
mensional quantity K4+„ is treated as a constant, the K4 
above also takes its currently measured value and is a 
constant rather than a variable) . 



The effective 4-dimensional curvature term is not 
canonical in Eq. (2); to make it so, let us take the con- 
formal transformation 

g,,u e^'^9t,i> (3) 
and choose the field (/) to satisfy 




Then we obtain the effective 4-dimensional gravitational 
action in the Einstein frame: 
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'^'Gravity — "2 ^ '^^"^ 
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(5) 



We shall make a further assumption that the extra di- 
mension(s) are homogeneous and isotropic, i.e. the met- 
ric of the extra space takes the following form: 



hab = diag(-6^ -6^ • • ■,-b'^), 
and then the action Eq. (5) could be rewritten as 
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by defining a new scalar field, the radion u: 



(T = 



K4 



log—. 
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(6) 



(7) 



(8) 



Now we turn to the matter sector of the effective 4- 
dimcnsional action, first considering the scalar fields. 
The action of a 4 + n dimensional scalar field (p (all the 
quantities with tilde are higher dimensional in this work) 
is given by: 



^i^+"X^/^G^ 



1 



(9) 



Since we are only considering the zero-mode theory and 
the zero-modes of the fields are independent of the extra 

dimensional coordinates (see Appendix A), we define a 
4-dimensional scalar field (which has the correct dimen- 
sion) (f from (fi: 



(10) 



ip{x) = ^/Voip{x,y). 



With this new field, the action Eq. (9) could be rewritten 
as 



d'xV\9\ 



Vo 



(11) 



where the new potential is related to the old one by 

U{^) = VoU{^). (12) 

Then the same conformal transformation Eqs. (3) and 
(4) transforms Eq. (11) into the following canonical form 
of the effective action: 



-J- 



\9\ i ^g^'^diJ.'^d^'f - exp 



— K 



n + 2 



(13) 



Note that from now on we will use k instead of Ki for 
simplicity. 

The same technique could be applied to gauge fields, 
whose higher dimensional action is given as 



'S'Ga 



= - j d'^xiTy^ 



2^ ^AB-: 



(14) 



where g is the A + n dimensional gauge coupling constant 
and are the corresponding gauge field strengths. 

Taking the following redefinitions of the zero-mode gauge 
field 



^1/ 



(15) 



and the conformal transformation Eqs. (3), (4), we finally 
obtain the effective 4-dimensional action in the Einstein 
frame as: 



auge 



■ j <^^a;v4sl^exp 



(16) 

One can obtain the effective action for the Dirac 
fermion field similarly: 



n + 2 



e^^^'zr'^Li^* + exp 



n + 2 



TO\1'\1' 



I 



(17) 



where Cj is the vierbein and F"^ are the Dirac matrices in the tangent space embedded in the 4-|-n dimensional flat 
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spacetime. The fermion mass m = y{H) is acquired from 
the Higgs mechanism with the Yukawa coupHng constant 
y and Higgs scalar field VEV (H) = v, which connect to 
their higher dimensional counterparts hy y — yj^fVa and 
V = ^/Vqv. 

To make the kinetic part of the fermion action Eq. (17) 
canonical, we rescale the field as 



and then by the conformality of the coupling of massless 
Weyl fermions, Eq. (17) becomes: 



il) = exp 



(18) 



— 



g\ \ e'^^tT'^D^^P + exp 



(19) 



r 



Our results above are equal to those of 37] when there 
are no large extra dimensions in their model. 



III. RADION DEPENDENCE OF 
FUNDAMENTAL CONSTANTS 

If there is a slow cosmological evolution of the extra di- 
mensional size between the time of BBN and now, then 
some or all of the fundamental constants in the parti- 
cle physics standard model will be changed and these 
changes may alter the results of the standard BBN. Be- 
cause the standard BBN is sensitively dependent on fun- 
damental constants, these changes, if exist, would be con- 
strained stringently by BBN. Furthermore, they have the 
potential of slightly modifying some aspects of the stan- 
dard BBN and improving the agreements between theo- 
retical calculation and observations. This section is de- 
voted to how the fundamental constants depend on the 
size of the extra dimensions (or equally the radion field). 



and in the next section we shall consider how BBN is 
influenced by these varying constants. 

As we are considering the system in the Einstein frame, 
the (4-dimensional) Planck mass will stay constant. Be- 
cause only dimensionless quantities such as the ratios of 
masses are physically significant, we shall take the Planck 
mass as a reference scale while expressing the variations 
of other quantities with dimension of mass such as Aqcd ■ 

The Higgs boson is a scalar field; therefore its radion 
dependence, as described by Eq. (13), is solely through 
an overall rescaling of the radion potential. As a result, 
the Higgs VEV, which is determined by minimizing the 
potential, will not be modified by the evolution of the 
radion. We shall take it to be a constant in the follow- 
ing calculation. Consequently, the radion dependence of 
fermion masses, as indicated by Eq. (19), should be due 
to the radion dependence of the 4-dimenioanl Yukawa 
coupling constant y. Let j/bbn and ?/now denote the 
Yukawa couplings at the time of BBN and now, and then 
from Eq. (19) we have: 
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(20) 



where we have used the definition p = Vbbn/^^^ow = 
^bbn/Vo- In deriving Eq. (20) we used the definition 
Eq. (8) of the radion field. The quantity p is useful for our 
purpose because it shows that the variations of Yukawa 
coupling constants (and other radion-relating quantities 
we shall consider later) depend only on the volume of the 
extra space, irrespective of how many extra dimensions 
there are. 

In the standard model, the Fermi constant Gp is not 
a real fundamental constant 13; rather, it could be ex- 



pressed as: 



9l 



(21) 



in which g2 is the coupling constant of the SU(2) gauge 
group and Mw is the mass of the weak gauge boson. In 
the extension of the standard model to UED scenario, 
relation Eq. (21) still holds, but with the Higgs field and 
weak gauge boson replaced by their corresponding zero 
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Kaluza-Klein modes (see Appendix A) and (72 by the 
effective 4-diniensional SU(2) gauge coupling constant. 
Therefore we conclude that the Fermi constant is unal- 
tered by the time evolution of the radion field. 

The effective gauge couplings could be read from 
Eq. (16) as: 
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Eq. (22) means that 



C^i.BBN -1 

— '■ = P , 

tti.NOW 



(22) 



(23) 



where i = 1, 2, 3 corresponds to the U(l), SU(2) and 
SU(3) gauge groups respectively, and the frequently used 
coupling constants at are related to gi by Ui — gf/An. 
The fine structure constant oem is obtained as a com- 
bination of the electroweak couplings, ~ So^T^/^ + 
a^^. Since ai and a2 have the same radion dependence 
described by Eq. (23), we have: 

OEM^ = p-\ (24) 
OEM, NOW 

Finally let us consider the influence of the radion evolu- 
tion on the strong coupling. The quantity more relevant 
to our calculation is the QCD scale parameter Aqcd, 
which is defined by the relation as(AQCD) — 00. Since 
c^s =0:3, the change of the strong coupling constant is 
given by Eq. (23). To obtain a relation between it and 
the change in Aqcd , we use the one-loop renormalization 
group equation for QCD governing the running of as: 



-11 



3"^ 



log^, (25) 



where E and Eq are any energy scales at which the values 
of ag are measured and np the number of quark flavors 
lighter than E. We shall choose the specified energy scale 
Eq to be the weak Z boson mass, il/z — 91.2 GeV . 
The present measured value of as at this energy scale is 
as(Mz) = 0.119 ± 0.002. Then the value of E at which 
ag^{E) 0, i.e. E — Aqcd, could be solved to be 



A 



QCD 



mbrric 



Ml 



exp 



27r 



9as(A/z 



(26) 



In deriving Eq. (26) we have made the assumption that 
Aqcd lies between the strange and charm quark masses 
TOs and TOc |38l l39l |. Note that unlike some previous 
works, we choose not to calculate the variations of the 
gauge couplings from the variation of a unified gauge 
coupling here, because we are concerned with the time- 
evolutions of the effective 4-dimenional gauge coupling 
constants at energy scales far below the inverse size of the 
extra dimensions, beyond which the field theory should 
be higher dimensional (or effectively the higher order KK 
modes should be taken into account). 

Using Eq. (23) and Eq. (26), we get the relation be- 
tween Aqcd.bbn and Aqcd,now: 



A 



QCD,BBN 



— p 27 exp 



A, 



QCD.NOW 



27r 



.9as,NOw(^^z) 



(27) 



IV. VARIATIONS OF QUANTITIES RELEVANT 
FOR BBN CALCULATION 



In this section we will briefly review the standard BBN 
theory (for more details see, e.g. 0j0,E3,E3|) ^"^^ "dis- 
cuss how its predictions are influenced by the variations 
of the fundamental constants considered in the previous 
section. Since some of these effects have been discussed 
in the existing literatures, we will not present them in 
details here. 

At very early times (T ^ 1 MeV or t <C 1 s), the en- 
ergy density of the universe is dominated by the photons, 
neutrinos and relativistic electron-positron plasma, with 
a negligible contribution from the baryons (mainly pro- 
tons and neutrons). These particles scatter frequently 
and are kept in thermal equilibrium. In addition, the 
rates of the weak interactions 



Ve + n ^ p + e'^ 
+ n ^ p + i7f, 

n ^ j5 + e^+z7g, (28) 

are far greater than the expansion rate of the universe, 
which is given by the Friedmann equation 



8nG 



"(?total) 



(29) 



where f?totai is the total energy density of the universe. As 
a result, there is also a chemical equilibrium among these 
particles so that the ratio between the number densities 
of neutrons and protons is 



exp 



T 



(30) 



in which TO„p = m„ — nip is the mass difference between 
neutron and proton. Eq. (30) is valid because in the 
standard BBN theory there is no lepton asymmetry and 
because of the charge neutrality of the universe |4Q]. 

As the temperature drops, the rates of weak interac- 
tions Eq. (28) decrease and become unable to keep neu- 
trons and protons in equilibrium. This begins to occur 
at a temperature of ~ 2 MeV [ij. After that, the actual 
n/p ratio will still be decreasing due to the free neutron 
decay and strong n—p reactions, until finally (essentially) 
all the neutrons have been processed into nulcei and their 
number density becomes constant. 

At least down to the temperature ~ 0.8 MeV the nu- 
clear reactions are capable of keeping the nuclei D, '^H, 
■^He, ^He in both kinetic and chemical equilibrium or the 
nuclear statistical equilibrium (NSE) with corresponding 
abundances |40l 
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YA=gA C(3) 



\A-1„ 



\(A-1) 



T 



(31) 



r 



where qa is the number of degrees of freedom of the nu- 
clear species (A, Z), tojv the nucleon mass, Ba the bind- 
ing energy of (A, Z), 77 the baryon-to-photon ratio and C, 
the Riemann zeta function. Because rj is of order 10^^", 
Eq. (31) says that the equihbrium abundances of the hght 
nuclei are very small at high temperatures. 

With the temperature decreasing further, the nuclear 
abundances provided by the reactions fall short of that 
required to maintain the equilibrium and depart from 
NSE, firstly for the heavier nuclei and later for lighter 
ones. The departure from NSE leads directly to the result 
that the back reaction rates become much smaller than 
the forward ones and are essentially switched off By 
T ^ 0.2 MeV, only the deuteron abundance is still held 
in NSE and the abundances for all other composite nu- 
clei are many orders of magnitude lower than their NSE 
values; some complex nuclei have been synthesized but 
the amounts are yet negligible. Then the heavier nuclei 
(''He, •^He, ^H) abundances begin to follow the deuteron 
NSE value. Finally, the deuteron bottleneck is passed 
at T ~ 0.1 MeV and a significant amount of deuterium 
is produced; the free neutrons are mostly rapidly assim- 
ilated into ^He, in which process some ^He and '^H are 
synthesized as the reaction ashes. These produced nuclei 
interact and lead to the production of tiny abundances of 
^Li, ^Li and ^Be, the last of which is finally turned into 
^Li by the electron capture process. 

At the end of BBN nearly all free neutrons are incor- 
porated in ''He, and so the final ''He abundance is very 
sensitive to the neutron-to-proton ratio at freeze-out and 
the duration between the freeze-out and the commence of 
BBN (for some neutrons decay in this period) . According 
to the review above, then, the final ''He abundance will 
be varied if at the time of BBN we have nonstandard val- 
ues for the neutron-proton mass difference m„p, the weak 
interaction rates and (or) the cosmic expansion rate, the 
last of which will be present in the case of a varying grav- 
itational constant or a varying energy density of some 
fiuid component in the universe. In contrast, the final 
yields of other complex elements depend strongly on the 
various relevant nuclear reaction rates. Thus we have to 
make clear how all these crucial physical quantities may 



J 



be modified with the varying fundamental constants as 
discussed in the above section. 



A. Neutron-proton Mass Difference 

The neutron-proton mass difference is given phe- 
nomenologically as [2^ 

runp = rUn - mp = rud - ruu - aEM^EM (32) 
with nid, rUu and asMWEM being respectively the down 
quark mass, up quark mass and electromagnetic self en- 
ergy difference, in which jtiem is determined by strong 
interactions and proportional to the QCD scale Aqcd- 
The quark masses are not exactly measured but could be 
calculated knowing the electromagnetic contribution as 
~ 0.76 MeV and the mass difference m„p ~ 1.293 MeV 
[4^ to be ~ 2.053 MeV. In the present model, since all the 
quantities appearing in Eq. (32) may change as the ra- 
dion field evolves and the variations of quark masses come 
completely from variations of the Yukawa couplings, the 
neutron-proton mass difference at the time of BBN is 
evaluated as 



"Ijip.BBN = 2.053 



2/BBN 
Z/NOW 



- 0.76 



ttEM.BBN Aqcd, BBN 
"EM, NOW Aqcd, NOW 



(33) 

which could be expressed in terms of p by using Eqs. (20), 
(24) and (27). 

The formulae Eqs. (32) and (33) are frequently 
used while considering rrinvS fundamental-constant- 
dependences, e.g. HllElllllE^- 



B. Weak Interaction Rates 

The weak rates of the neutron- proton inter-conversions 
can be estimated using the Fermi theory. The n — > p and 
p — > n rates are the summation of the rates of three 
forward and backward reactions in Eq. (28) respectively 
and could be expressed as [ioj : 



r(n ->p) 



A 



de 



A 



de 



e(e-(z)2(e2_ 1)1/2 



[1 + exp(-eze)] {1 + exp [(e - q)z^]} 
e(e + Q)2(e2_ 1)1/2 



[1 + exp(eze)] {1 + exp [-(e -I- q)z^]}' 



(34) 
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T{p n) 



A 



de- 



e(e-g)2(e2_ 1)1/2 



A 



[1 + exp(eze)] {1 + exp [{q - e)z^]} 
I [1 + exp(-eze)]{l+ exp [(e + 5)^1/]}' 



de- 



(35) 



where we have used 



'Tip 



dimensionless quantities q — 
~ m^/T^ and = 'nie/T 
T respec- 



with rrie being the electron mass, and 
tively the temperatures of the neutrinos and the elec- 
tromagnetic plasma. A is a normalization factor de- 
termined by the requirement that at zero temperature 
T{n p + e~ + Vg) = where t„ is the neutron life- 
time. This means that: 



A- 



'A(g)- 



^ cx Ghnl 



where 



de e{e-qy{e^-l 



a/2 



(36) 



(37) 



We have shown above that the Fermi constant Gp is 
independent of the radion field, and so the variations of 
weak rates at BBN from their current values originate 
from the changes of nie and q, or equivalently from the 
changes of the Yukawa coupling, the fine structure con- 
stant and the QCD scale Aqcd- It is then straightfor- 
ward to calculate the neutron lifetime and weak interac- 
tion rates at different values of p with the aid of Eqs. (20), 



(24), (27), (33) and (36). In our calculation we have done 
the explicit numerical integrations of Eqs. (34) and (35) 
in the BBN code to take into account the effects of an 
evolving radion. For this purpose A is obtained by firstly 
computing Anow using Eq. (37) and the measured value 
of Tn^NOWj a-iid then calculating Abbn with the aid of 
Eq. (36). We have not evaluated explicitly the various 
corrections to the weak rates which are left for fu- 
ture work; instead we use the Wagoner's approximation, 
i.e., reducing the n ^ p rates by 2%. The difference 
between these two treatments is estimated by Lopez and 
Turner to be ^ 0.5% for ''He, which lies well within 
its observational uncertainty (see below). 



C. Expansion Rate of the Universe 

Because we are working in the Einstein frame, the 
gravitational constant will not change. However, the 
change of the electron (and also positron) mass does 
induce a variation of the energy density stored in the 
electron-positron plasma because the energy density and 
pressure both depend on the electron mass |45j |: 



J 



Qe+ + Qe 



9 °° 1 



3 1 

-Ksinze) + -Ki{nze) 



Pe+ +Pe- = —^t X! 



(_)"+! 1 

cosh{n(f>e ) K2 {nze ) , 



(38) 
(39) 



in which 0e is the electron chemical potential and Ki (i 
= 1, 2, 3) are the hyperbolic Bessel functions '4d|. In the 
actual calculations we truncate the infinite summations 
in Eqs . (38) and (39) to retain only the first 5 terms as 
in [3q . 

Since this energy density is an important part in the 
total energy budget of the universe, we expect that the 
expansion rate will change also according to Eq. (29), 
though this effect is not significant. Note that the vari- 
ation of Aqcd will cause the masses of all the nucleons 
and nuclei to be different from their present-day {i.e. the 
standard BBN) values, acting as another source of non- 
standard energy density. However, because the contribu- 
tion of baryons to the total energy density is negligible, 
here we simply ignore this effect. 



D. Nuclear Reaction Rates 



The nuclear reactions relevant to BBN could be di- 
vided roughly into the charged-particle-induced and the 
neutron-induced ones, the former of which are charac- 
terized by the Coulomb barrier penetration since the two 
charged particles must overcome their Coulomb repulsion 
to approach each other and interact. In the authors 
analyze the possible infiuences of a varying fine structure 
constant on the cross sections of these charged-particle 
reactions. By writing the cross section as 



aiE) . ^exp 



— 27raEM^i^2W ^ 



(40) 
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in which S{E) is the astrophysical factor, a slowly vary- 
ing function of energy, Zi {i =1, 2) are the charges of the 
reacting nuclei and ji is the reduced mass of the system, 
they assume that all the dependence on oem lies in the 
exponential term in Eq. (40). Then various thermonu- 
clear reaction rates are calculated using the fitting for- 
mulae of Smith, Kawano & Malaney ^41] (hereafter SKM) 
and expressed in terms of the variation of asM- Later, 
NoUett & Lopez Q improve the work of 0] by includ- 
ing several effects not considered there. In the present 
work we use the original proposal of and include the 
improvements by in evaluating the reaction rates. 

In addition to the influences from q;em, the variation 
of the reduced mass ^ in the exponent of Eq. (40) could 
also change the cross section. This is because, with the 
same energy, a different particle mass implies a differ- 
ent particle velocity, and, since the penetration depends 
on the velocity, it consequently implies a different pos- 
sibility of penetrating the Coulomb barrier. Meanwhile, 
these cross sections may be changed with the modified 
nuclear binding energies and a few reactions have reso- 
nance terms which are considered to be dependent on the 
strong couplings. However, lacking explicit expressions 
for the coupling dependence of these effects, we choose 
to neglect them as in 26j and treat the main contribu- 
tion as coming from the exponent in Eq. (40); an esti- 
mation of the errors by neglecting these effects will be 
presented in Sec. VI (it is reassuring to observe the fact 
that, since nuclear masses can be well approximated as 
proportional to Aqcd and so /i (x Aqcd, has a simi- 
lar varying trend as that of Q!em and in fact several or 
many times larger in magnitude. Thus the exponent in 
Eq. (40) may be rather significant). The thermonuclear 
rates taking into account the Aqcd dependence are not 
difficult to obtain in forms similar to those in us- 
ing the thermo-averaging method of Fowler, Caughlan & 
Zimmerman |47| . 

There may also be dependences of the binding energies 
of the complex nuclei on the strong coupling, which are 
not understood very well but will also modify the inverse 
reaction rates. However, we have checked that these ef- 
fects are typically rather small and negligible (at least 
for our allowed range of parameters given below in Sec- 
tion VI) except for the reaction p{n, D)7, which will be 
discussed below. 

In contrast to the charged-particle-induced reactions, 
for the neutron-induced reactions there are no Coulomb 
penetrations, and their strong coupling dependences 
should be considered in the context of the theoretical 
calculations of the cross sections. Consider the reaction 
p(n, D)7, the first step of BBN. The cross section of this 
reaction is poorly determined experimentally j43 |. For- 
tunately we have a rather good theoretical understand- 
ing for it and whenever a comparison is possible, the 
measurements show excellent agreements with theoreti- 
cal calculations. In this work, we adopt the calculation 
with the low energy effective field theory without pions 
[isL l49| . The authors compute and express the p(n, 0)7 



cross section as a function of the deuteron binding energy 
Bd, the nucleon mass Af^r, the scattering length in the 
singlet channel and so on, concentrating on the energy 
range relevant for BBN. The rate we obtain by thermo- 
averaging this theoretical cross section |4S| is virtually 
identical to the data-fitting results of plL l50l | . To esti- 
mate the fundamental constant dependence of the rate 
we, following ,2S | ■ assume the parameters which appear 
in this theoretical formula and have dimensions [Length]" 
to be proportional to rn~"'. The deuteron binding en- 
ergy Bd is important because it appears not only in the 
cross section formula but also in the rate of the reverse 
reaction 7(D, n)p (and in fact it is the latter that is dom- 
inant). Recent studies of the pion-mass-dependence of 
the nuclear potential show a strong variation of Bd with 
the pion mass m^r jSll l52j |: particularly, in contrast to 
previous results, these authors show that Bd is in fact 
a rfecreasmo rather than increasing function of m,r- As 
noticed in U H^, although the calculations in |^ Is^ l 
have large uncertainties, we could safely approximate the 
deuteron binding energy Bd by 



Bd 



,BBN 



= B, 



d.NOW 



(r + 1)- 



'71,1, BBN 



TO7r,NOW 



g(TO7r,BBN) 



(41) 



within the narrow range we are interested in, i. e. a small 
variation of m^. In Eq. (41) the parameter r varies from 
6 to 10 according to the calculations of [sills^j . For the 
present work we choose r = 10 js^l- We have checked 
the case of r = 6 and found that there is only a ~ 0.05% 
increase in the best fitting value of p while the features 
discussed below are exactly the same. 

To relate Eq. (41) to the fundamental constants dis- 
cussed in the previous section, we use the Gell-Mann- 
Oakes-Renner relation [s^ for pion mass: 



f^ml = (m„ + m,d){qq) cx to. 



1/2 



(42) 



in which the coupling of the pion to the axial current /^r 
and the quark condensate (qq) are proportional to Aqcd 
and Aq(-,p) respectively and iriq denotes quark mass. Note 
that in |52l| the authors, while quoting their results in 
terms of the pion mass m^, were varying the ratio of the 
quark mass to Aqcd , and so they indeed showed a range 
of Bd vs TO^ which is effectively a function 8] (using 
Eq. (42)) as 



Bd,BBN — Aqcd,BBN5' 



TOq,BBN 



A 



QCD,BBN 



(43) 



where the function g' is related with function g in Eq. (41) 
through a change of variable. In Ref. 8] Eq. (43) is used 
to obtain the relation between Bd and TO„p, both of which 
depend on the quantity Aqcd; here we adopt the same 
relation to estimate the variation of Bd in the present 
model. There is also a weak dependence of Bd on the 
fine structure constant asM because of the small electro- 
magnetic contribution (0.018 MeV) to Bd H^, and we 
also have included this effect in our calculation. 
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FIG. 1: (color online) The primordial abundances of ^He, D and ^Li as a function of the baryon-to-photon ratio 77, for various 
values of p as indicated by the legend in the lower right panel. Also shown are the observational allowed ranges of these 
abundances (the horizontal cyan lines) and the WMAP- indicated range of rj (the vertical blue lines). 



With the help of Eq. (43) and the theoretical cross 
section for p(ri, D)7, it is then straightforward to obtain 
the radion-dependence (p-dependence) of the thermo- 
averaged rate for this reaction. 

For the other two neutron-induced reactions 
'^Re{n,p)^li and ''Be(n,p)'^Li, there is no theoreti- 
cal calculation like that for p{n, D)7, and so we simply 
use the data-fitting results and neglect their coupling 
dependences. The reaction ^He(ri,,p)^H influences the 
•^He output which is not used to constrain the parame- 
ters, and so the ignorance of its coupling dependences is 
non-essential. But the ''Li abundance which we use to 
reduce the parameter space does depend on the reaction 
^Be(n,p)^Li and this reaction exhibits a broad resonance 
which may shift if the fundamental couplings vary; in 
Sec. VI we shall give a simple estimation on the error 
due to neglecting its coupling dependences. 



V. NUMERICAL RESULTS 

We have incorporated all the above-discussed effects 
into the standard BBN code by Kawano [s^. For the 
nuclear reactions, we re- fit the updated rates given in the 



NACRE compilation [s^ (and a later work of Cyburt, 
Field & Olive js^) in the forms as in SKM to compute 
the impacts of varying asM and fi. Then we compare 
this result with that obtained by simply using the SKM- 
fitted rates and find that they are essentially the same. 
Keeping this in mind, in the following discussions we will 
only use the results derived from the SKM rates. 

In FIG. 1 we plot the abundances of the light nuclei 
"'He, D and ^Li including the effects of an evolving radion 
field, which is characterized by the ratio between the ex- 
tra space volumes at the time of BBN and at present, 
p. It is apparent that if the size of the extra dimensions 
is larger at the BBN era {p > 1), then there will be an 
increase in the deuterium output and decreases in the 
"^He and ^Li (for ^Li we only consider the larger-7y-case in 
consistent with the WMAP result J,,]) yields compared 
with the standard BBN. The behavior of the ''He out- 
put is mainly the consequence of two effects originated 
respectively from the weak interactions and the nuclear 
reaction p{n, D)j. On one hand, a p- value larger than 
1 implies a smaller uem at the time of BBN, leading to 
a larger rn„p (from Eq. (33)) and smaller neutron den- 
sity at freeze-out, thus finally to a smaller ^He output 
(note that as discussed above, the dependence on the re- 
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FIG. 2: (color online) Similar to FIG.l, but here we neglect the effects coming from the reaction p{n, D)^ (mainly the deuteron 
binding energy B^)- See the text for explanations. 



duced mass jj has the same trend and strengthens this 
effect). Meanwhile, the weak interaction rates will be 
larger than those in the standard case, again decreasing 
the final "'He abundance. On the other hand, the larger 
p-value will cause a smaller deuteron binding energy Bd 
(from Eq. (43)) than in the standard BBN, and so the nu- 
cleosynthesis will commence at a later time and become 
less efficient, producing less '*He and ''Li while leaving 
more D unprocessed (there is an extra decrease in the 
forward rate of p{n, D)"f which is also due mainly to B^, 
but we have checked that its influence is small compared 
with the effect of the later commencing BBN). Both of 
these two effects work in the same direction for ^He so 
that the final ^He yield is very sensitive to p. As for D and 
^Li, the effects of a larger p other than that of Bd come 
mainly from the Coulomb barrier: the charged-particle- 
induced reaction rates will be less suppressed since both 
Q!EM and p are smaller in this case, consequently more D 
is processed and more ^Li is produced. To see this point 
more explicitly, we plot in FIG. 2 the BBN yields if all 
the effects discussed in the previous section except that 
of the p{n, D)j reaction are included. FIG. 2 agrees with 
the trend when only the effects of changing qem are con- 
sidered, such as that in [2>^3- However, the inclusion of 
varying Bd (as in FIG. 1) changes these abundances dra- 



matically, indicating the important role of the deuteron 
binding energy in BBN 0,11, [l3|. 

Our results as shown in FIG. 1 suggest that the dis- 
crepancy between the standard BBN theory and WMAP 
observations tends to be reduced if p is greater than 1 
Pt. Using the WMAP-imphed value of 0, the stan- 
dard BBN can reproduce the observed D abundance, but 
not the ''He and ^Li abundances; an rj smaller than the 
WMAP value is needed to obtain the observed "'He and 
(or) ^Li abundances. While this discrepancy may result 
from some systematic errors in the ^He and ''Li observa- 
tions, the possibility of nonstandard BBN or new physics 
cannot be excluded. In fact, there have been many at- 
tempts of solving this discrepancy from a theoretical as- 
pect, for example, by nonstandard expansion rate [s^ . 
varying fine structure constant and lepton asymme- 
try etc. These effects separately cannot solve the 
whole inconsistency and a combination of them is needed, 
6.17. |17| | (there are also works which use a single quan- 
tity to account for the inconsistency, such as a varying 
deuteron binding energy and the Brans-Dicke cos- 
mology with a varying A term |6ll|). But from a theoret- 
ical viewpoint, such combinations may not be necessary 
since, as mentioned in the introduction, the variation of 
one fundamental constant will often be accompanied by 
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FIG. 3: (color online) The 68% (solid) and 95% (dashed) contours from the observational data on D (upper left), *He (upper 
right) and ^Li (lower left). The lower right panel shows the 68% contours plotted together to make the comparison explicit. 



the variations of some other constants; these changes to- 
gether would play the role of several effects combined. 
Because the co-variation of the fundamental constants is 
rather model-dependent, a general study is not practical. 
Our present work serves as a specific example and shows 
how the variations of several fundamental constants are 
correlated with one parameter p and driven by the same 
physics - in this case, the variation of the volume of extra- 
dimensional space. 

In FIG. 1 the horizontal lines show the la ranges of 
the measured ^He, D and ''Li abundances. We use the 
observational results of Olive, Steigman and Walker ma] 
for '^He and the most recent work of Kirkman et al. |63| 
for D: 



T^obs 



0.238 ±0.005, 
p/jj)obs ^ 2.78t|]-j^ X 10"^ 



(44) 
(45) 



For Lithium we use the value given in the recent work of 
Bonifacio et al. 1641. 



('Li/H) 



obs 



2.19t°:^^ X 10 



-10 



(46) 



which is larger than previous results but still a factor 
of 2 smaller than the standard BBN -f WMAP result. 



The vertical lines in FIG. 1 represent the Icr range of r] 
allowed by the CMB analysis [l|: 
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?/WMAP = (6.14 ±0.25) X 10" 



(47) 



In order to see to what extent the quantity p can differ 
from 1 , we next carry out the likelihood analysis and find 
the confidence contours in the two-dimensional parame- 
ter space (r],p). For this purpose, we adopt the semi- 
analytical method introduced in [63 and then slightly 
generalized in • The likelihood function is simply 

given as L oc exp [— x^/2] in which the chi-squared is 
calculated by: 



i,3 



with Wij the inverse covariance matrix: 



2,th 



2. obs r 



(48) 



(49) 



2,th 



and Yl {Y° ) are the theoretical (observational) abun- 
dances of the i-ih. element. In Eq. (49) dj^j 
matrix calculated by 



is the error 



2,th 



\T.['^^in)-Y^{V-,)] [Y,{Tt)-Y,{T, 



(50) 
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FIG. 4: Joint constraint on the parameter space from D+^He+'^Li; the gray and light gray regions denote the 68% and 95% C. 
L. contours respectively. The point with the maximum likelihood, which lies around (5.9 x 10"^", 1.0011) on the (??, p) plane, 
is also shown (the white circle). The vertical lines represent the la range of r) determined by WMAP and the dot-dashed 
horizontal line is the case of the standard BBN (p = 1). 



in which F^, is the fc-th nuclear reaction rate plus (mi- 
nus) its ItT uncertainty; the summation is over all the 
most relevant reactions. When i = j, Eq. (50) simply 
gives the theoretical uncertainty of the i-th nuclear abun- 
dance. 

In FIG. 3 we show the 68% and 95% C. L. contours 
on the rj — p plane calculated using the modified Kawano 
code. Although neither D nor ^He alone could constrain 
the parameters considerably, we could see from the lower 
right panel of this figure that the combination of them 
reduces the allowed parameter space significantly, a sim- 
ilar conclusion as that in the previous works 0, [2a | . The 
WMAP-implied value of -q is seen to lie in the allowed 
parameter space of anyone of D & ''He, D & ^Li and 
''He & ^Li. We have also plotted the joint constraint 



on the parameters from D-l-^He-l-'^Li in FIG. 4 with the 
WMAP-implied la range for 77, Eq. (47), indicated by 
the vertical lines. It can be seen that our best fitting 
value (the white circle) lies just on the edge of the la 
range of ?7wmap while the whole la WMAP range falls 
into our 68% C. L. contour. On the other hand, the stan- 
dard BBN {p ~ 1) cannot be consistent with the WMAP 
result and is disfavored at the 68% C. L. though it still 
lies inside our 95% C. L. contour. The allowed range 
of p and ?7 by our calculation is 0.9999 ^ p ^ 1.0021, 
4.94 X 10-1° < V ^ 7.05 x 10"^° at the 68% C. L. and 
0.9992 < p < 1.0027, 4.46 x lO'^" < 77 < 7.99 x 10'^° at 
the 95% C. L., of which the ranges of p could be trans- 
formed into the allowed variations of the fine structure 
constant in the present model as: 



- 2.1 X 10-=^(-2.7 X 10-=^) < "EM.BBN-aEM.NOW < ^ ^q_3(o g ^ ^Q-3y (5^) 

Q!EM,NOW 

I 



This constraint is looser than that obtained in [26j but 
more stringent than that in 0]. 



VI. DISCUSSIONS AND CONCLUSIONS 



In summary, we have considered in this work the im- 
plications of a slowly evolving radion field on the fun- 
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damental constants and the Big Bang Nucleosynthesis 
predictions in the context of the universal extra dimen- 
sions model. We have included in our modified BBN 
code various effects of the varying constants, of which 
the most important ones being the neutron-proton mass 
difference, the deuteron binding energy, the free neu- 
tron lifetime and the nuclear reaction rates. These ef- 
fects may act in either the same or the opposite di- 
rections and lead to distinct and complex behaviors of 
the light nuclei yields. We have also calculated the chi- 
squared as a function of p, the ratio between the ex- 
tra space volumes at the time of BBN and now, and 
r], the baryon-to-photon ratio; this analysis then gives a 
bound on the variation of the fine structure constant as 
-2.1 X 10-3 < AaEM/ttEM < 0.1 X 10-3 at 68% C. L. and 
-2.7x 10-3 < AaEM/aEM ^ 0.8x IQ-^ at 95% C. L. The 
allowed range for AaEivi/aEM in our result is much larger 
than that obtained from quasar absorption systems [20j| . 
which is AckEM/aEM = (-0.72 ±0.18) x 10-^. However, 
because BBN probes physics up to redshifts of order 10^° 
while the quasar absorption systems are at z = 1 ~ 3, 
there should be no apparent inconsistency. 



It is suggested that a larger extra dimension size at 
the time of BBN may help reduce the discrepancy be- 
tween the BBN theory and observation, and here we 
make several comments about this conclusion. Firstly, 
the expression for the neutron-proton mass difference 
rrinp (Eq. (32)) is phcnomcnological but reasonable, even 
though it may not be definitive |^ we here use it 
because it is the most-commonly-used and we find no 
competing ones. Secondly, there is an uncertainty for 
the parameter r in the deuteron binding energy formula 
(Eq. (41)) derived by the effective field theory, but as 
we stated above, changing r from 10 to 6 only induces a 
small variation of the numerical results while having no 
influence on the qualitative features we have discussed. 
Thirdly, our discussion on the variations of nuclear reac- 
tion rates follows that in [2^ and we have neglected sev- 
eral effects due to the lack of theoretical understandings 
for them. One example is the influence of a varying Bd on 
the reaction rates of D(p, 7)3He, D(D,n)3He, D(D,p)^H, 
D(3H,n)''He, D(3He,p)*He, which are dominant or sub- 
dominant in the destruction of D. Lacking an explicit 
expression for this effect, we choose not to include it in 
the modified code. Instead we can make an estimation 
about the error introduced by neglecting it, by assum- 
ing the cross sections of these reactions to scale with the 
deuteron radius and thus Bd, i.e., a oc 1/ {misiBd) 
Adopting this parametrization in the code, we find that 
the outputs of D, ^He and ''Li are changed by — ^ 9%, 
< 0.12% and < 5% respectively in our interested range of 
(77, p) (see FIG. 4), which lie well within the correspond- 
ing la observational uncertainties (^ 16% for D, ~ 2% 
for ^He and ~ 21% for ''Li). Another example is the 
influence of resonances in some key cross sections such 
as those of D(3h, n)4He, D(3He,p)4He and ''Be(n,p)''Li. 



These resonant terms are generally of the form [42 



(a|t;|)„s = F(r)exp 



E 



(52) 



in which -F'(r) is a function of the temperature and E 
is the resonance energy. Without a reliable calculation 
of how these terms depend on the strong coupling, we 
roughly estimate the corrections from them by assuming 
E cx Aqcd- Then using the fitting values of E in these 
cross sections the modified code we find that the 

D and ""Li abundances decrease by < 0.14% and < 1.4% 
in the same range of (ry, p) as above, while the '^He out- 
put stays essentially unchanged. These values again fall 
well in the corresponding observational uncertainties. Al- 
though the points discussed in this paragraph would not 
change our results much, they do refiect the important 
role these effects play in reducing possible errors, and in 
any case the precisions of the BBN constraints will im- 
prove with future progress in particle and nuclear physics 
relevant for these topics. 

It is also interesting to look at whether the allowed 
range for p in our constraint would imply additional ef- 
fects that are too large to be dangerous for other ob- 
servations, such as the existence of bound di-proton/di- 
neutron and the stability of ^He. It is well known that 
the bound di-proton will open a rapid channel for the 
hydrogen fusion 67] and thus be catastrophic to the star 
lifetimes. The binding condition for di-nucleon systems 
was considered by, 6.17. Dent and Fairbairn |3^ . and a 
rough criterion for these systems to be stable was a de- 
crease of Aqcd to ~ 1/10 of its present day value (see 
[3^ for more details), which is much larger than our con- 
strained range for the gauge coupling variations. The 
stability of ^He is also important in the consideration of 
varying fundamental constants as, if stable, ^He would 
fill the mass-5 gap in the BBN nuclear reaction chain, 
drastically enhancing the production of ''Li. This issue 
was discussed by, e.g. Flambaum and Shuryak Q, and it 
was found there that the limit from ^He binding cannot 
compete with that from B^. For example, even the more 
stringent criterion to avoid a stable '^He quoted in Ref. |^, 
namely = S{mq/ Aqcd)/ {mq/ Aqcd) ^ -0.05, is weh 
satisfied by our 2a range of p. 

As a final point, we comment on the choices of ob- 
servational ''He and ''Li abundances. The value given 
in Eq. (44) is taken from Ref. and is certainly not 
the only observational result of Yp. Other data analy- 
sis give different resuhs, such as Y^^"" = 0.2421 ± 0.0021 
in Ref. 68] and Y^^"" = 0.249 ± 0.009 in Ref. 69]. It is 
apparent from FIGs. 1 and 3 that, if these larger abun- 
dances are used, the '^He contour will shift rightward. 
Meanwhile, as stated above, the ''Li abundance as taken 
from Ref. |6J] is significantly larger than previous esti- 
mates such as CLi/H)"''" = 1.23 ± 0.06tj^ (95% C. L.) 
given by Ryan et al. [t^- From FIGs. i and 3 we see 
that, with the latter result adopted, the right branch of 
the ''Li contour (which we are interested in) will shift 
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leftward. Consequently the inconsistency between BBN 
and WMAP as discussed above would be less reduced. 
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APPENDIX A: A FIVE DIMENSIONAL TOY 
MODEL OF UED 

In this appendix we review briefly the main ingredients 
of the minimal universal extra dimension model relevant 
for our discussion. In particular we show that the UED 



model reduces to 4D SM at low temperatures. For sim- 
plicity we will assume only one extra dimension in this 
toy model, which is assumed to be flat and compactified 
on an orbifold Si/Z^. The fifth dimension is character- 
ized by its coordinate X4 = y G [0, 2ttR], where R is the 
radius of the fifth dimension (Note that here the conven- 
tion is different from that in Section II; see for more 
details) . 

As a minimal generalization of SM, the 5-dimensional 
model contains three generations of quarks and leptons, 
denoted by Qi, Ui, Di, Li and Ei respectively where 
i = 1, 2, 3 is the generation index (For example, the zero 
modes oi Qi, Ui, Di, Li and Ei correspond respectively 
to the (u,d)L, Ufi, da, (i/e,e~)i, and in the 4D SM). 
Then the 5-dimensional Lagrangian in the matter sector 
could be expressed as: 

-^5, total = Lq + Lp + L}i + Ly, (Al) 

in which the subscripts G, F, H and Y represent the 
gauge, fermions, Higgs and Yukawa sectors of the model 
respectively and are described as: 



Lg = 

Lf = 

La = 
Ly = 



- V— Tr 



F^^''{X)FrAB{X) 



X V^Da 



Q{X),U{X),D{X),L{X),E{X) 



DaH{X) D^H{X) - fl^H+{X)H{X) -~X H+{X)H{X) , 
YuQ{X)U{X) ■ i(j2H*{X) -f YdQ{X)D{X)H{X)YeL{X)E{X)H{X) + H.c. 

I 



(A2) 

(A3) 

(A4) 
(A5) 



In the above F^^ denotes the 5-dimensional gauge field 
strength associated with the SU(3)c x SU(2)w x U(1)y 
group and Da the corresponding covariant derivative; 
F^'s are the 5-dimensional anti-commuting gamma ma- 
trices defined as Ta = {7^,7^^} with 7^ {^i ^ 0,1,2,3) 
being the 4D Dirac matrices and 7'* = —i^^j^j'^j^ (for 
the definition of gamma matrices in more than 5 dimen- 
sions see J33I)- As in the standard model the fermion 
masses are generated by the Higgs mechanism, the 5D 
Higgs doublet H{X) has the same form as in 4D, ji and A 
being its potential parameters; the Yukawa couplings K's 
are 3x3 matrices with mass dimension —1/2. Notice that 
in Eq. (A3) there is an implicit summation over the 3 gen- 
erations of quarks and leptons and that the capital Latin 
letters A and B run from to 4: X^ = , ,y. 



In the 4-dimensional effective theory the fields could be 
written as an infinite summation of Kaluza Klein (KK) 
modes and the zero KK modes are believed to play the 
role of the SM particles at low energies. However, the 
KK expansion itself does not necessarily lead to an exact 
equality of the zero KK modes to SM; to obtain such an 
equality the KK decomposition should satisfy that, (1), 
only left- (right-) handed component of each weak dou- 
blet (singlet) is even under the orbifold projection and 
(2), the redundant zero modes of the 5th components of 
gauge fields, A4 , should be eliminated. The KK expan- 
sions with appropriate boundary conditions which satisfy 
the above requirements are given as: 



HiX) 



V2ttR 



H^-°\x) + %/2y ij(")(x) 



1 



X COS' 



ny 
R 



(A6) 
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MX) 
Q{X) 

U{X) 

b{x) 
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The factor y/2 above is due to the different normahza- 
tions of zero and higher order modes in the KK tower 
and will disappear if we run the summation over both 
positive and negative values of the KK numbers n. It's 
apparent from the above decompositions that the zero 
modes are independent of the extra dimension coordinate 
y, a fact which is expected because the standard model 
should be purely 4-dimensional. The higher KK modes, 
however, generally depend on y; and they will acquire 
additional masses of order n/_R, which will be accessible 
only with higher enough energies because of the small- 
ness of R (l/R should be larger than several hundreds 
of GeV 113 )■ Thus up to the energy range relevant for 
BBN (0(1) MeV) the world acts as 4-dimensional and 
the SM description shall be safe. 

Lastly let us turn to the Higgs sector of this model. 
The fermion masses are only generated after the spon- 
taneous symmetry breaking (SSB); after that, the Higgs 
field acquires a vacuum expectation value which is deter- 
mined by minimizing the Higgs potential 



in which /i^ < so that the potential has nontrivial min- 
ima [iJl ■ In the 4D effective theory the mass of the n-th 
Higgs KK mode is given by m^^^''^ = /i^ + 11^ / R^ . There- 
fore, if the size of the extra dimensions is small enough 
such that |/i| < 1/i?, then m^^'^ for n > 1 is positive, 
and only the neutral zero KK mode obtains a nonzero 
expectation value [tJ- In this condition the zero Higgs 
doublet KK mode could be expanded as: 



75 



(A15) 



(A14) 



where /i'-^-' is the physical Higgs zero mode 0^°)+ and x'"^ 
are the zero modes of Goldstone bosons; v is the (zero- 
mode) Higgs expectation value. At low energies as in our 
case, the effects of the higher order Higgs KK modes are 
negligible and will not be considered here. 
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